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Abstract 

We propose a method to determine the isoscalar KN scattering length on the lattice. Our method 
represents the generalization of Liischer's approach in the presence of inelastic channels (complex 
scattering length). In addition, the proposed approach allows one to find the position of the S- 
matrix pole corresponding the the A(1405) resonance. 
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1. The antikaon-nucleon scattering amplitude is of fundamental importance in nuclear, particle and 
astrophysics, see e.g. Refs. [1-4]. In particular, the KN system at threshold provides an interesting 
testing ground of the chiral dynamics of QCD with strange quarks due to the A(1405) resonance just 
below the scattering threshold. In fact, experimental information on the K~p scattering length from 
scattering data and kaonic hydrogen level shifts is contradictory, as first stressed in [5] and further 
elaborated on in Refs. [6-8]. A clarification is expected from the upcomimg SIDDHARTA experiment 
at DA$NE, that intends to remeasure kaonic hydrogen with unprecedented accuracy and is expected 
to give further constraints to the isoscalar and isovector kaon-nucleon scattering lengths from the first 
measurement of the energy spectrum of kaonic deuterium. 

On the theoretical side, effective field theory methods in various disguises are employed to pin down the 
KN scattering length. The scattering information is usually analyzed in terms of unitarized versions 
of chiral perturbation theory, that lead to a dynamic generation of various resonances, in particular 
the A(1405). Such schemes have been worked out over the years by various groups at various levels of 
sophistication (see, e.g. [6-9]). While an impressive amount of data (cross sections, threshold ratios, 
mass distributions, etc.) is described in such approaches with good precisison, unitarization of course 
introduces some unwanted model-dependence. On the other hand, the extraction of the KN scattering 
length from kaonic hydrogen is firmly rooted in non-relativistic bound-state effective field theory and 
thus is devoid of the above-mentioned model-dependence (for a recent review, see [10]). However, it 
entirely rests on the availablity of precise kaonic atom data. Unfortunately, the existing data from 
DEAR [11] and KEK [12] are conflicting. It would therefore be most welcome to have another tool at 
hand that would allow one to determine this fundamental quantity. 

As we will argue in this letter, lattice QCD provides such a framework. As first shown by Liischer, 

finite volume simualations of the energy levels of two-particle states can give access to scattering 
information [13, 14]. The idea of Liischer is very elegant and simple. For two-particle states that are 
well separated from bound states or resonances in the given channel, the 1/L expansion of the energy 
levels takes the generic form (here, L is the size of the box with volume L x L x L) 



where a is the desired scattering length and ci,C2 are pure numbers (see below). The method has 
e.g. been used to extract the tttt, ttK and KN S-wsve scattering lengths from lattice data [15]. Note 
that an alternative proposal to extract the scattering length from the two-particle wave function is 
e.g. given in Ref. [16]. 

However, for the extraction of the KN scattering length, a generalization of this scheme is called for 
since there is a strong channel coupling between KN and Svr, the latter channel having its threshold 
about 100 MeV below the opening of the KN one. In addition, the appearance of the A(1405) just 
between these two thresholds further complicates the picture. All these features can be captured by a 
two-channel Lippmann-Schwinger equation. As we will show in the following, a suitable formulation of 
this equation in the finite volume allows for an unambigous extraction of the complex- valued isoscalar 
K~N scattering length. Note also that our method is similar to the approach adopted in Ref. [17], 
where the problem was treated within the potential scattering theory. In this paper, we use instead the 
language of the non-relativistic effective field theory (NR EFT), which enables one to systematically 
address the effects of particle creation/annihillation and relativistic corrections. 

2. To set the stage, we consider a two-channel Lippmann-Schwinger (LS) equation in NR EFT in the 
infinite volume. Note that we are using a covariant version of the NR EFT, considered in Refs. [18]. 
The channel number 1 refers to KN and 2 to Stt with total isospin 1 = 0. The resonance A(1405) is 
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Figure 1: Complex s-plane with the Svr and cuts along the real axis and the location of the 
A(1405) resonance. 

located between two thresholds, on the second Riemann sheet, close to the real axis (these thresholds 



are defined by st = {rriN + Mk)"^ and s[ = (ms + Mtt)^) W% Consider first energies above ^A'^ 
threshold, s > {mj\f + Mk)"^- The coupled-channel LS equation for the T-matrix elements Tij{s) in 
the dimensionally regularized NR EFT reads (we only consider 5-waves here) 

Til = -f^ii + Hii iqiTu + Hu iq2T2i , 

T21 = H21 + H21 iqiTu + H22 iq2T2i , (2) 

with qi = A^/2(s,m^,M|.)/(2^), q2 = \^/'^{s,m^,Ml)/{2y/s) and \{x,y,z) stands for the Kallen 
function. Furthermore, the Hij{s) denote the driving potential in the corresponding channel. Con- 
tinuation of the center-of-mass momentum qi below threshold (ms + M^^)'^ < s < {rriN + Mk)"^ is 
obtained via (see Fig. [1] for the corresponding analytical structure) 

zgi - -Ki (3) 

The resonance corresponds to a pole on the second Riemann sheet in the complex s-plane, its position 
can de determined from the secular equation 

A(s) = l + K^Hii- 4 H22 - Kf4 {H11H22 - Hl^) (4) 

with Kf = -(-A(sH,m2^,Mi))i/2/(2^) and = (-A(s^,, m|, M2) V2)/(2^). The energy and 
width of the resonance are then given by y/sR = Eft — iTji/2. The KN scattering length is related to 
the amplitude Tn at s = st = {rriN + M^^)"^ via 

-rr ( \ u ( N , i92(st)(iii2(st))^ ... 
aii = Tii[st) = Hii{st) + — TTTTT^ • ^ 

1 - n2{St) H22{St) 

Thus, to pin down its complex value, we need to determine the three real quantities Hn, H12, H22 at 
s = St appearing in Eq. ([5]) 

3. We now consider the same problem in a finite volume. The rotational symmetry is broken to a 
cubic symmetry so that the infinite volume version of the LS equation Eq. ([2j) takes the form (we 
consider only S'-waves here, neglecting the small mixing to higher partial waves. The mixing can be 
easily included at latter stage, see e.g., Ref. [19].), 

2 2 

Tu = Hii Zoo(l; A;i) -ffiiTii Zoo{l; H12T21 , 



^2\ = -^21 — 1=^ -^00(1; ^1) H21TU — 1=^ -^00 (1; ^2) -^^22221 



(6) 



**In the following, we work in the isospin limit and thus do not resolve the further splitting of these thresholds. 
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Here, we have neglected the terms that vanish exponentially at a large L. The secular equation that 
determines the spectrum can be brought into the form 
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This is rewritten as 



Zoo{l;H)F{s,L)=0 , 
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(8) 
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-4>{k2) + nvr , n = 0, 1,2, 



arctan 



Zoo(l;A:i) ' 



(9) 



with 



tan 5{s, L) = q2{s) F{s, L) . (10) 

5{s,L) is called the pseudophase. It is a function of the energy ^/s and the level index n, 5n(s) = 
5(s,L„(s)). 

The dependence of the pseudophase on s and L (or, equivalently, on the level index n) is very different 
from that of the usual scattering phase. Namely, the elastic phase extracted from the lattice data 
by using Liischer's formula is independent of the volume modulo terms that exponentially vanish at 
a large L. Further, the energies where the phase passes through tt/2 lie close to the real resonance 
locations. In contrast with this, the pseudophase contains terms which are only power suppressed at 
a large L. Moreover, it contains the tower of resonances which are not related to the dynamics of 
the system in the infinite volume and merely reflect the existence of the discrete energy levels in the 
"shielded" channel. 

Measuring the pseudophase on the lattice can be used to determine the KN scattering length. It can 
be directly seen from the expression of the pseudophase, which depends on real functions Hn, Hi2,H22- 
Extracting these from the data, we then find the scattering length by using Eq. ([5]). Note that in the 
expression for the scattering length we need Hij[s) evaluated at threshold s = st- We shall however 
demonstrate below that replacing Hij{s) by Hij[st) in certain observables, related to the pseudophase, 
introduces very small correction, since the effective range term proportional to {s — st) is suppressed 
by as compared to the leading order result. To be specific, we consider the following three 
observables: 

1. For some chosen value of n, we measure the value of the pseudophase (5(st; L{st)) = 5t at threshold 
St = {rriN + Mk)"^ and Et = y/sl (see Fig. [3]). On the other hand, we may express 5t through 
St in the following way. At threshold, Zqo is singular. 



H'l j at 5 



Zoo(l;A;2 



1 1 



4tt k"^ 



ill) 



4 



so that 

F{s,L)U^s, = H22{st)-HUst)/Hii{st) (12) 
iB.n5{st]L{st)) = q2ist) {H22{st)-Hf2ist)/Hu{st)) =q2ist)I{st) . (13) 

Thus, measuring 6t, we may extract the combination H22 — -^^12/^11- 

2. Suppose that tan5{s; L{s)) is infinite at s = = and L = = L{s3) (see Fig. [2] for a 
specific representation of the pseudophase based on a two-channel K-matrix model described 
below). This occurs at the energy where the denominator of Eq. ([8]) vanishes 

l-^^oo(l;fc?(s3))i^ii(53) =0. (14) 

We solve this equation by expanding both iJii(s) and ^00 (1; ^i('S)) in Taylor series in the vicinity 
of s = St 

Hu{s) = Huist) + ql{s) H[^{st) + 0{q^) (15) 



and 



^ Zoo{l;k') = ^^ + '^ + ^{cl-C2) + Oik') , (16) 



with 

1 ^r^ 1 



limV -4^ = -2.837297... , C2 = c? - ^ V 4r = 6.375183 ... . (17) 
Substituting Eqs. and (flUj) into Eq. ([HD, we obtain 



qAs3) = 13— (^l + -i-^ + -2^2- + 0(l3jJ • (18) 

This means that measuring the value of s, where the tan 5(s;L(s)) becomes infinite, we may 
extract Hn at s = st- Note that the effective-range term, which contains H[i{st), contributes 
first at ©(L-S). 

3. Similarly, suppose that tan5{s,L) = at s = S2 = -E'i L = L2 = L{s2) (see Fig. [2]). Using 
the same technique as just described, we obtain: 

where 

G{st) = Hnist) - Hl2{st)/H22{st) • (20) 
Therefore, measuring S2; we get Hu — Hf2/H22- 

Using finally Eq. ([5]), we can express the scattering length in terms of the three quantities Hu, I 
and G, all taken at s = sj 

au = Hu[st)-\ : — / X r/ n tj / n n • (21) 

1 - iq2{st) I{st)Hii{st)/G{st) 

This is the central result of this letter. 
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Figure 2: The pseudophase 6. The energy Ei at which the pseudophase passes through n/2 corre- 
sponds to the A(1405). 



Finahy, note that in the analysis of the lattice data it may be more convenient to directly fit the 
explicit expression of the pseudophase given in Eq. ([8]) to the measured values on the lattice around 
s = St, replacing Hij{s) by Hij{st) and considering Hij{st) {ij = 11,12,22) as three independent 
fitting parameters. From the above discussion one may expect that such a fit will lead to the precise 
determination of these parameters. The effective range terms can be neglected since their contribution 
is suppressed by three powers of L. In this case, one does not need to measure the pseudophase in the 
whole interval between S2 and S3. 

4. Given the parameters Hij determined from fitting to the pseudophase, the position of the pole on 
the second Riemann sheet of the complex variable s, which corresponds to the A(1405)-resonance, can 
be determined from the secular equation dH). We expect that replacing Hij{s) by Hij{st) allows one 
to locate the pole position at a reasonable accuracy. 

5. In order to demonstrate the above-described proposal in practice, we have investigated a coupled- 
channel model with an explicit A(1405) resonance located at Rey^jR = 1406 MeV and — 2Imy^SR = 
50MeV. Effective range terms are neglected. The matrix elements Hij are taken equal to 

Hu = -1.47573 fm , Hu = 0.91581 fm , H22 = -0.34159 fm . (22) 

This corresponds to an = aQ^K^N) = (—1.26 + i0.70)fm. The resulting first four energy levels as 
a function of M„L are shown in Fig. [3l We notice that the lowest level (n = 1) does only show a 
moderate volume dependence in the interval considered, quite in contrast to the excited ones with 
n > 2. For M^-L ~ 2 ... 3 the ground state level flattens around E = 1406 MeV that corresponds to 
the A(1405). It is clear that, for this reason, the lowest level can not be used for the extraction of the 
KN scattering length. The excited levels show a more complicated behavior in this interval of L. At 
the first glance, these levels exhibit the so-called avoided level crossing somewhere between 1430 MeV 
and 1440 MeV. In the elastic case, such a behavior of the energy levels signalizes the presence of a 
narrow resonance near this energy. However, this is not the case here. The peculiar behavior of the 
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Figure 3: Energy levels for the two-channel model with an explicit A(1405) resonance in the finite 
volume. The avoided level crossing which is observed at the energies between 1430 MeV and 1440 MeV 
is not related to the physical resonance in the infinite volume but reflects the presence of the KN 
threshold. For comparison, we plot the energy levels levels for the non-interacting two-particle systems 
ttS (dashed lines) and KN (dotted lines). 

excited energy levels is caused by the opening of the KN threshold. At higher energies, the picture 
repeats ~ an avoided level crossing emerges, if the KN system has a discrete eigenvalue at this energy 
in a finite volume. If the volume changes, the avoided level crossing moves (in difference to the avoided 
level crossing corresponding to the "true" resonance). For L — > cxd the bifurcation lines accumulate at 
threshold s = st- In this limit, the scattering amplitude is not analytic at s = (unitary cusp). 

In Fig. [2] gives the pseudophase derived from the second (n = 2) energy level. It shows the expected 
behavior. First, it crosses 7r/2 at ^/s = Ei, very close to the mass of the A(1405). Then, it passes 
vr at ^/s = y/s2 = E2, close to the threshold where its value is (5j > vr. At £'2, the tangent of the 
pseudophase vanishes at since qf{s2) < 0, we can conclude that G{st) > 0, cf. Eq. (fT9|) . Finally, the 
value of 37r/2 is reached at ^/s = ^/s^ = E^. Here, (?i(s3) > and consequently Hii{st) < 0. This can 
be deduced from Eq. (fT9]) after the substitutions S2 — > S3 and L2 — > L3. 

6. In this letter, we have generalized Liischer's algorithm for the extraction of the scattering length 
from the finite-volume energy spectrum measured on the lattice. The modified algorithm applies to 
the case when the scattering length is complex due to the presence of the open channel(s) below 
threshold. In the case of the KN scattering with total isospin / = 0, the scattering length can be 
determined by measuring the volume dependence of the first excited level around the threshold energy. 
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